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ABSTRACT 

We estimate the power spectrum of HI intensity fluctuations for a sample of 8 
galaxies (7 dwarf and one spiral). The power spectrum can be fitted to a power 
law Pm{U) — AU" for 6 of these galaxies, indicating turbulence is operational. 
The estimated best fit value for the slope ranges from ^ —1.5 (AND IV, 

NGC 628, UGC 4459 and GR 8) to 2.6 (DDO 210 and NGC 3741). We 

interpret this bi-modality as being due to having effectively 2D turbulence 
on length scales much larger than the scale height of the galaxy disk and 
3D otherwise. This allows us to use the estimated slope to set bounds on the 
scale heights of the face-on galaxies in our sample. We also find that the power 
law slope remains constant as we increase the channel thickness for all these 
galaxies, suggesting that the fluctuations in HI intensity are due to density 
fluctuations and not velocity fluctuations, or that the slope of the velocity 
structure function is 0. Finally, for the four galaxies with "2D turbulence" 
we find that the slope a correlates with the star formation rate per unit area, 
with larger star formation rates leading to steeper power laws. Given our small 
sample size this result needs to be confirmed with a larger sample. 
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1 INTRODUCTION 



Evidence has been mounting in recent years that turbulence plays an important role in the 
physics of the ISM as well as in governing star formation. It is believed that turbulence is 
responsible for ge nerating the hierarchy of s tructures present across a range of spatial scales 
in the ISM (e.g. lElmegren fc Scalo 



2004a 



Elmegren fc Scald l2004bl ) . In such models the 



ISM has a fractal structure and the power spectrum of intensity fluctuations is a power law, 
indicating that there is no preferred "cloud" size. 

On the observational front, power spectrum analysis of HI intensity fl uctuations is an im- 



porta nt technique to probe the structure of the neutral ISM in galaxies (ILazarian 



1995 



Oev 



2OO2I ). The power spectra of the HI i ntensity fluctuations in our own galaxy , the LMC and the 



SMC all show power law behavior ( 



2000 



Elmegreen et al 



2001 



lent medium. Similarly, 



Crovisier &: Dickev 



Stanimirovic et al 



Westpfahl et al.i 



1983 



Green 



1993 



Deshpande et al 



1999D which is a characteristic of a turbu- 



19991) showed that the HI distribution in several 



Krakow et al 



(11982) have estimated the 



galaxies in the M81 group has a fractal distribution, 
power spectrum of optical intensity fluctuation for the galaxies M 81, M 51 and NGC 1365 
using digitized images in differ ent bands. Th e powe r spectra were found to show both long 



and short range order. Further, 



Willett et al. 



(I2OO5I ) used the Fourier transform power spec- 



tra of the V and Ha images of a sample of irregular galaxies to show that the power spectra 
in optical and Ha pass-bands are also well fit by power laws, indicating that t here is no 



characteristic mass or luminosity scale for OB associations and star complexes. iRoy et al. 



( 120091 ) have estimated the power spectrum of the supernovae remnants Cas A and Crab. The 
power spectrum is found to be power law, indicating the presence of magneto-hydrodynamic 
turbulence. 



Recently iBegum et al.l (120061 ) have presented a visibility based formalism for determining 



the power spectrum of HI intensity fluctuations in galaxies with extremely weak emission. 



This formalism is very similar to that used for analyzing rad io- interferometric o 



the Cosmic Microwave Background Radiation anis otropies (IHobson et al. 



Dservations of 



1995 



White et al. 



1999 



Hobson fc Maisinger 



2002 



Dickinson 



20041 ). This formalism was first applied to a 



dwarf galaxy, DDO 210. Interestingly, the HI power spectrum of this extremely faint, largely 
quiescent galaxy was found to be a po wer law with the same slope (—2.8 ± 0.4) as that 



observed in much brighter galaxies. We (iDutta et al 



20081 ) have applied the same technique 



to an external spiral galaxy NGC 628 (M 74) and found a different slope (—1.7 ±0.2). In this 
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Galaxies 


DDO 210 


NGC 628 


NGC 3741 


UGC 4459 


GR 8 


AND IV 


KK 230 


KDG 52 


(1) Type 


Dwarf 


Spiral 


Dwarf 


Dwarf 


Dwarf 


Dwarf 


Dwarf 


Dwarf 


(2) Distance (Mpc) 


1.0 


8.0 


3.0 


3.56 


2.1 


6.7 


1.9 


3.55 


(3) Log[SFR] (Moyr-l) 


>-5.42 


-0.1 


-2.47 


-2.04 


-2.46 


-3.0 


>-5.53 


>-5.1 


(4) Angular extent 


5' X 4.5' 


13' X 12' 


27' X 12' 


4.5' X 4' 


4.2' X 4' 


15' X 8.5' 


3' X 2' 


3.5' X 3' 


(5) o-o{,3(km s-l) 


6.5(1.0) 


8.5(1.5) 


- 8.0 


9.0(1.6) 


9.0(0.8) 


~ 8.0 


7.5(0.5) 


9.0(1.0) 


(6) iHl 


26.0° 


6.5° 


68.0° 


30.0° 


27.0° 


55.0° 


50.0° 


23.0° 


(7) References 


3, 6 


1,8 


4, 6 


5,6 


2, 6 


7, 6 


5, 6 


5, 6 



T able 1. 1-lKa.rnpliuis fc Briggj lll992l) 2-lBeeum fc Chengalurl fcOOSf) , 3- lBeeum fc Chengaluj ll2004h. 4-lBegurn et al.l l l2007^ . 
5- iBegum et al. I ll2006h . 6- lBegum et all ||2008^ ■ 7- Chengalur et al. 2009 (in preparation')8 - iPersic fc Rephaelil (120071) 



paper we use the same formalism to measure the power spectrum of HI intensity fluctuation 
for a sample of 6 more dwarf galaxies. For completeness, we also include the existing results 
for DDO 210 and NGC 628, and search for correlations with other galaxy properties, e.g, the 
star formation rate (SFR). Any detected correlation would provide insights into the nature 
of turbulence in the ISM. In section 2 we discuss the galaxies in our sample. In section 3. we 
discuss the power spectrum estimator used in our work. In section 4, we perform simulations 
to support the analytical results obtained in section 3. The method of analysis is discussed 
in section 5. In the last section we present our results and conclusions. 



2 THE GALAXY SAMPLE 

We briefly discuss some of the relevant properties of the galaxies in our sample. A few 
of the galaxy parameters are also summarized in Table [1] which contains - Row(l): type 
of the galaxy; Row(2): distance in Mpc; Row(3): Log[SFR] in MQyr~^ measured from Ha 
emission; Row(4): angular extent of the galaxy calculated from Moment map at HI column 
density 10^^ atoms cm~^; Row(5): velocity dispersion (km s~^), Row(6): inclination in HI 
and Row (7): references. 



DDO 210 



DDO 210 is the faintest {Mb ~ —10.9) relatively close (at a distance 950 ±50 kpc; 



Lee et al. 



19991 ) gas-rich member of the Local Group. The HI disc of the galaxy is nearly face-on. On 



large scales, the HI distribution is not axisymmetric; the integrated HI column density 
contours are elongated towards the east and south. No Ha emission was detected indicating 
a lack of on-go ing star formatio n in the galaxy. Distance to this galaxy is estimated to be 



950±50 kpc in 



Lee et al. 



fll999h . 
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NGC 628 

NGC 628 (M74: ) is a nearly face-on SAfs Ic spiral galaxy with an inclination angle in the 
range 6° to 13° JKamphuis fc BriggJ[l992) . It has a yery l arge H I disk extending out to more 



than 3 times the Holmberg diameter. 



Elmegreen et al. 



(120061 ) have found a scale-free size 



and luminosity distribution of star forming regions in this galaxy, indicating turbulence to 
be functional here. The d istance to this galaxy is uncertain with previous estimates ranging 



from 6.5 Mpc to lOMpc. 



Sharina et al. 



(119961 ) estimated a distance of 7.8 ± 0.9 Mpc from 



the brightest blue star in the gala xy. This distance e stima te matches with an independent 



photo metric distance estimate by 



Sohn fc Davidge I (119961 ). In a recent study 



Vinko et al 



(120041 ) inferred the distance to be 6.7 ± 4.5 Mpc by applying the expanding photosphere 
method to the hyper novae SN2002ap. In this paper, we adopt the photometric distance of 
8.0 Mpc for this galaxy. 



NGC 3741 

NGC 3741 is a nearby dwarf irregul ar galaxy (Mr ^ —13.13) with a gas disk that extends 



to ~ 8.8 times the Holmberg radius (IBegum et al. 



2005 



20071 ). The galaxy is fairly edge-on 



with kinematical inclination varying from ~ 58° — 70°. NGC 3741 appears to have a HI bar 
and is very dark matter dominated with a dark to luminous mass ratio of ~ 149. Further, 
this galaxy is undergoing significant star formation in the center. An i nterplay between th e 



neutral ISM and star formation in this galaxy is studied in det ail in 



The rotation curve flattens beyond 300" to a value ~ 50 kms ^ ( Gentile et al. 



Begum et al 



(l2007h 



2007). They 



Karachentsev et al 



find th at the ISM in this galaxy shows radial motion of ~ 5 — 13 km s~^. 
(I2OO4J ) have estimated the distance to this galaxy as 3.0 ± 0.3 Mpc using the tip of the red 
giant branch (TRGB) method. 



UGC 4459 

The faint dwarf {Mb ~ —13.37) galaxy UGC 4459 is a member of the M 81 group of galaxies. 
It is fairly isolated from its nearest neighbor UGC 4483 at a projected distance of 3.6° (~ 223 
kpc) and a velocity differe nce 135 km s~^. UGC 4459 is a relatively metal poor galaxy, with 



12 + log(0/H) ~ 7.62 (iKunth fc Osthn 



2OOOI ). The optical appearance of this galaxy is 



dominated by bright blue clumps, which emits copious amount of Ha, indicating high star 
formation. The velocity field of UGC 4459 shows a large scale gradient across the galaxy 
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with an average of ~ 4.5 km kpc~^, though this gradient is not consistent with that 
expected from a systemat ic rotating disk. This galaxy has a TRGB distance of 3.56 Mpc 
(iKarachentsey et aLll2004j ). 



GR 8 



'his is a faint (Mr ^ —1 2.1) dwarf irregular galaxy with very unusual HI kinematics 
(iBegum fc Chengalur ( 20031 )). The HI distribution in the galaxy is very clumpy, and shows 



substantial diffuse, extended gas. The high density HI clumps in the galaxy are associated 
with optical knots. In optical it has a patchy appearance w ith the emis sion dominated by 
bright blue knots which are sites of active star formation (IHodgd 119671 ). Both radial and 
circular motions are present in this galaxy. It also possess es a faint extended emissio n in Ha. 
The distance to this galaxy is estimated to be 2.10 Mpc (IKarachentsev et al.ll2004j ). 



AND IV 



This is a dwarf irregular galaxy with a rnodera te surface brightness (/ly ~ 24) and a very 
blue colour (V — I < 0.6) (Ferguson et al.ll2000l) . It is at a projected distance of 40' from the 



center of M 31 (IMerrett et al 



20061 ) and has a very low ongoing SFR of ~ O.OOlMQyr ^. 



Its HI disk extends to ~ 6 times the Holmberg diameter and also shows large scale, purely 
gaseous spiral arms Chengalur et al. 2008 (in prepa ration). The distance to this galaxy is 



estimated to be 6.7 ± 1.5 Mpc (IFerguson et al. 



2000). 



KK 230 



KK 230, the faintest(MB ~ —9.55) dwarf irregular galaxy in our sam ple, lies at the periphery 
of the Canes Venatici I cloud of galaxies (IKarachentsev et al.ll2004l ). The velocity field shows 
a gradient in the direction roughly perpendicular to th e HI and optical major axis with 
a magnitude of ~ 6 km s^^ kpc~^ (IBegum et al.ll2006l ). There is no measurable ongoing 



star formation in this ga. 



Karachentsev et al 



axy as inferred from the absence of any detectable Ha emission. 



(120041 ) have found a tidal index of —1.0 indicating the galaxy to be fairly 



isolated. They have estimated the TRGB distance to be 1.9 Mpc. 
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W(U) 




Figure 1. W{U) for top-hat and exponential window functions for = 1'. The vertical hne marks 9^ ^ . Note that the FWHM 
of the two window functions are nearly the same and both window functions have a very small value for U > 9^^ . 



KDG 52 

KDG 52 (also called M 81DwA), another member of the M81 group, is a faint dwarf galaxy 
{Mb ~ —11.49) with a clum py HI distribution in a broken ring surrounding the optical 



emission (IBegum et al. 



20061 ). The HI hole is not exactly centered around the optical emis- 



sion. This galaxy does not have any d etectable ongoing star for mation. The distance to 
this galaxy is estimated to be 3.55 Mpc (IKarachentsev et al.ll2004l ) derived from the TRGB 
method. 



3 A VISIBILITY BASED POWER SPECTRUM ESTIMATOR 

The specific intensity of HI emission from a galaxy can be modeled as 

iu{e) = wM [h + 6h{e)\ . (1) 

Here 9 is the angle on the sky measured in radians from the center of the galaxy. We assume 
that the galaxy subtends a small angle so that 9 may be treated as a two dimensional 
(2D) planar vector on the sky. The HI specific intensity is modeled as the sum of a smooth 
component and a fluctuating component. Typically, Iv{9) is maximum at the center and 
declines with increasing 9. We model this through a window function Wy{9) which is defined 
so that Vrjy(O) = 1 at the center and has values 1 > Wy{9) > elsewhere. This multiplied 
by lu gives the smooth component of the specific intensity. For a face-on galaxy, the window 
function Wi,{9) corresponds to the galaxy's radial profile. 

Consider a galaxy of angular radius 9q. In our analysis we consider two different models 
for the window function of such a galaxy. In the "top-hat" model it is assumed that the 
specific intensity has a constant value within a circular disk of radius 9o, and it abruptly 
falls to outside. The window function, in this Heaviside step function 
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W,{9) = Q{9o - 9) (2) 
In the "exponential" model the window function has the form 



WM = exp [-^j ■ (3) 

where it falls exponentially away from the center. 

It is also possible to use Wu{9) to define a normalized window function Wi^{9) such 
that / (P9W^{9) = 1. The second moment of 9 defined as / (P9 9^W^ {9) provides a good 
estimate of the angular extent of the window function. Here we have used the condition that 
this should have the same value 9q/2 for all models of the window function, to determine 
the width of the exponential window function in terms of 9q. Throughout we have assumed 
^ 1 radian. 

We express the fluctuating component of the specific intensity as Wv{9) 5Iv{9). Here 
5Iv{9) is a stochastic fluctuation which is assumed to be statistically homogeneous and 
isotropic. The HI emission traces these fluctuations modulated by the window function 
which quantifies the large-scale HI distribution. In this paper we would like to use radio- 
interferometric observations to quantify the statistical properties of the fluctuations 5Ip{9). 
These fluctuations are believed to be the outcome of turbulence in the inter-stellar medium. 

The visibility Vu{U) recorded in radio- interferometric observations is the Fourier trans- 
form of the product of the antenna primary beam pattern Ay{9) and the specific intensity 
distribution of the galaxy. 

V,{U) = Jd9e~''^'^-'A,{9)h{9) (4) 

Here U refers to a baseline, the antenna separation measured in units of the observing wave- 
length A. It is common practice to express U in units of kilo wavelength (kA). Throughout 
this paper we follow the usual radio inteferometric convention and express U in units of kilo 
wavelengths (kA). 

The angular extent of the galaxies that we consider here is much smaller than the primary 
beam, and the effect of the primary beam may be ignored. We then have 

V.(f7) = W{U% + W{U) ® 6X{U) + K{U) (5) 

where the tilde ~ denotes the Fourier transform of the corresponding quantity and (S> denotes a 
convolution. In addition to the signal, each visibility also contains a system noise contribution 
MviU) which we have introduced in eq. ([5]). The noise in each visibility is a Gaussian random 
variable and the noise in the visibilities at two different baselines U and U is uncorrelated. 
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The Fourier transform of the normahzed window functions are 
^ Ji{27reoU) 



and 



1 



(6) 



(7) 



[1 + 7rXUy3f^ 

for the top-hat and exponential models respectively. Here Ji(x) is the Bessel function of 
order 1. These functions, shown in Figured], have the property that they peak around U = 
and fall off rapidly for U ^ 9q^. This is a generic property of the window function, not 
restricted to just these two models. At baselines U ^ 9q^ we may safely neglect the first 
term in eq. ([5]) whereby 

V.(f7) = W{U) ® + K{U) (8) 

We use the power spectrum of HI intensity fluctuations Pm{U) defined as 

{6X{U)6V 0')) = 610 - U') Pm{U) (9) 

to quantify the statistical properties of the intensity fluctuations. Here, 5|)(f/ — U') is a two 
dimensional Dirac delta function. The angular brackets here denote an ensemble average over 
different realizations of the stochastic fluctuation. In reality, it is not possible to evaluate 
this ensemble average because a galaxy presents us with only a single realization. In practice 
we evaluate an angular average over different directions of U . This is expected to provide 
an estimate of the ensemble average for a statistically isotropic fluctuation. 
The square of the visibilities can, in principle, be used to estimate PmiU) 



and this has 



(10) 



Lazarian 



1983 



Green 



1993 



W^U) ® Phi{U) + {\ K{u) I') 

Deen used in several earlier studies (ICrovisier fc Dickey 
19951 ). This technique has limited utility to observations where the HI signal in 
each visibility exceeds the noise. This is because the last term (| My{U) p), which is the 
noise variance introduces a positive bias in estimating the power spectrum. The noise bias 
can be orders of magnitude larger than the power spectrum for the faint external galaxies 
considered here. In principle, the noise bias may be separately estimated using line-free 
channels and subtracted. In practice this is extremely difficult owing to uncertainties in 
the bandpass response which restrict the noise statistics to be estimated with the required 
accuracy. Attempts to subtract out the noise-bias have shown that it i s not possible to do 
this at the level of accuracy required to detect the HI power spectrum (lBegumll2006l ). 
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10 r, ^ ^ ^ , 



0.1 



P( U, AU) 



0.001 



le-05 I ' ' ' 1 

1 10 

U (kX) 

Figure 2. This shows the convolution (eqn.[TTJ of a power law Phi(C^) = U~'^'^ with the product of two exponential window 
functions (eqn. [7]l with 9o = 1' ((ttSq)^^ ~ IkA). The bold solid curve shows the original power law and the thin solid curve 
shows Piii{U , AU) for AU = 0. The dashed curves are for| AU |= 0.2, 1.0 and S.OkA respectively from top to bottom. The 
direction of AU is parallel to that of U. Note that PmiU, 0.2) Piii{U, 0) and the correlation between two different visibilities 
is small for | At? |> {7r6»o)"^ 

The problem of noise bias can be avoided by correlating visibilities at two different 
baselines for which the noise is expected to be uncorrelated. We define the power spectrum 
estimator 

Pm{U,AU) = ( V.(f7)V:(f7 + Af7) ) 

= J <fU' W^{U - U') W*{U + AU- U') Pm{U') (11) 

Since W{U) falls off rapidly for U > 9^-^ (Figure [I]), the window functions Wu{U — U') 
and W*{U + Af/ - U') in eqn. ([11]) have a substantial overlap only if |Af/| < (7r6'o)"^. 
Visibilities at two different baselines will be correlated only if |A[7| < (tt^^q)^^) and not 
beyond (Figure [2]). In our analysis we restrict the difference in baselines to |Af/| ^ {-kOq)^^ 
so that Wi,{U + Af/ — U')'^Wu{U — U') and the estimator PHi(f^, Af/) no longer depends 
on Af/ (Figure [2]). We then use the visibility correlation estimator 

PHi(f/) = ( V.(f7)v:(f7 + A?7) ) 

= j d^U' I W,{U - U') |2 PHi(f/') . (12) 

The measured visibility correlation pHi(f^) will, in general, be complex. The real part is 
the power spectrum of HI intensity fluctuations convolved with the square of the window 
function. A further simplification is possible at large baselines U ^ ^, provided | Wy{U) p 
decays much faster than the variations in PHi(f/). We then have 

¥ni{U) = C Pm{U) (13) 

where C = / | Wu{U) p d^U is a constant. 

We use the real part of the estimator PHi(f/) to estimate the power spectrum PHi(f^)- 
Our interpretation is restricted to the U range U ^ (tt^^q)"^ where the convolution in eq. 
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-3 -2.5 -2 -1.5 -1 -0.5 

a 



Figure 3. This shows how Um changes with a for an exponential window function with 6o = 1'. Note that the discontinuities 
seen in the plot appear to be genuine features and not numerical artifacts, though the cause of these features is not clear at 
present. 
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Figure 4. Effect of window function modifies the power spectrum differently for different power law exponents. 

( TT2I) does not affect the shape of the power spectrum and eq. ( |T3|) is a vahd approximation. 
In order to estimate the U range where this approximation is vahd we have numerically 
evaluated eq. (fTTj) assuming PmiU) to be a power law PYniJJ) = AU°^. Figure H] shows 
the results for both the top-hat and the exponential models with 6q = 1' (9q^ = 3.4 kA) 
and a = —1.5 and —2.5 which roughly spans the range of slopes that we encounter in our 
galaxy sample. Using Um to denote the value (in kA) where the deviation from the original 
power law is 10%, we find that for the top-hat and exponential models respectively Um has 
values (3.1,5.3) for a = —1.5 and (4.8,13.0) for a = —2.5. Note that Um depends on two 
parameters, namely a and 6q. Figure [3] shows how Um changes with a for the exponential 
model with 9o = 1'. For other values of 9q we scale the value of Um in Figure [3] using 
Um oc 9q^. For a given power law index, the estimator Phi(^) gives a direct estimate of the 
power spectrum for U > Um- 

The estimator Phi(?7) also has a small imaginary part that arises from the HI power 
spectrum because the assumption that Wu{U + AU — U) ^ W„{U — U) is not strictly vahd. 
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U (k X) 



P(k) 



le+07 



le+05 



le+03 



10 




100 



Figure 5. The 3D and 2D power spectrum P{k). The simulated HI power spectrum PHl(f^) for the thin (A) and thick (B) 
disk without the radial profile are also shown. The P{k) and k values (left and bottom axes) have been arbitrarily scaled to 
match the Phi(C^) and U axes (top and right). The power-law fits are shown by solid lines. The different curves have been 
plotted with arbitrary offsets to make them distinguishable. 



We use the requirement that the imaginary part of PHi(f^) should be small compared to the 
real part as a self-consistency check to determine the range of validity of our formalism. 

The real and imaginary parts of the measured value of the estimator PHi(f^) both have 
uncertainties arising from (1.) the sample variance and (2.) the system noise. The sample 
variance is [PniiU) / y/T^]"^ , where A''^; is the number of independent estimates. We assume 
that in the u-v plane the HI signal is uncorrelated beyond 1.5 times the FWHM of |iyj(f/)p. 
We use this to determine Ne from the u-v coverage of our observations. The system noise 
variance is a^/Np where is the variance of the individual visibilities in the data and Np 
is the number of visibility pairs that contribute to each visibility correlation. We add both 
these contributions to determine the 1 — a error-bars. The reader is referred to Section 3.2 
of 



AhetaL 



(120081) for further details of the error estimation. 



4 SIMULATION 

We have carried out simulations to asses the impact of the overall galaxy structure on our 
estimates of the power spectrum of HI intensity fluctuations. The starting point is a three 
dimensional (3D) 512^ cubic mesh, labeled using coordinates r, on which we generate a 
statistically homogeneous and isotropic Gaussian random field h{r) with power spectrum 
P{k) = Ck"* . We have used 7 = —2.5 throughout our simulations, the results can be easily 
generalized to other 7 values. The 3D power spectrum is shown in Figure [5l For reference 
we also show the two dimensional (2D) power spectrum of h{r) evaluated on a 512^ planar 
section of the cubic mesh. As expected, the 2D power spectrum also is a power law with 
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Figure 6. The simulated HI power spectrum for the face-on thick disk with (B) no radial profile, and (C),(D) with radial 
profile using R = 0.35,0.7 respectively. (E) is same as (C) with the disk tilted at 60°. The different curves have been plotted 
with arbitrary offsets to make them distinguishable. 
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Figure 7. The simulated HI power spectrum for the face-on thin disk with (A) no radial profile, and (F),(G) with radial profile 
using R = 0.35,0.7 respectively. (H) is same as (F) with the disk tilted at 60°. The different curves have been plotted with 
arbitrary offsets to make them distinguishable. 

P{k) (X k'^'^^. In all cases we have generated five independent realizations of the Gaussian 
random field, and averaged the power spectrum over these to reduce the statistical uncer- 
tainties. 

In our simulations we use the Gaussian random field h{r) as a model for the 3D HI 
density fluctuations that would arise from homogeneous and isotropic 3D turbulence. The 
h{r) values on a 2D section through the cube serves as a model for the HI density fluctuations 
in the limiting situation where we ignore the thickness of the galaxy and treat it as a 2D disk. 
Note that the resulting h{r) is a statistically homogeneous and isotropic Gaussian random 
field on the 2D section. We use this to represent the density fiuctuations that would arise 
from homogeneous and isotropic 2D turbulence. 

We next embed a galaxy in the middle of the 3D cube. The overall, large-scale structure 
of the galaxy is introduced through a function G(r) so that the HI density at any position is 
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G{f)[hQ + h{f)]. Here G{f)ho is the smoothly varying component of the galaxy's HI density 
and G{f)h{r) is its fluctuating component. It is assumed that the observer's line of sight is 
along the z axis. The HI density G{r)[hQ + h{r)] is projected on the x — y plane. We interpret 
the projected values as the HI specific intensity lu{0) in the plane of the sky (eq. [1]). The 
512^ mesh in the x — y plane of our simulation corresponds to 4' x 4' on the sky. 

We first consider a face-on galaxy, and use G{r) = exp(— z^/^;^). This incorporates only 
the finite thickness of the disk which is characterized by the scale-height Zh and it ignores 
the galaxy's radial profile in the plane of the disk. Note that, for z < Zh this function 
closely matches the function sec hj—z^/zl) which is also used to model the scale-height 
profile (IBinney fc Merrifieldlll998l ). For the scale- height, we have used the values Zh = 4: and 
32 mesh units which we refer to as the "thin' and "thick" disk respectively. Figure [5] shows 
the power spectrum of the simulated HI specific intensity lui^) for both these cases. We 
find that the HI power spectrum of the thin disk has a slope —2.5, same as the 3D power 
spectrum, at large U. There is a break at f/ ~ 35 kA which corresponds to l/irzh, and the 
slope is ~ —1.9 at smaller U. We interpret this change in the slope in terms of a transition 
from 3D fluctuations on scales smaller than the disk thickness (large U) to 2D fluctuations 
on scales larger than disk thickness (small U). We expect the slope to approach —1.5, the 
2D slope, at very small U. Note that in our simulation it is not possible to evaluate the 
slope at very small U values where the sample variance is rather large. We find that the HI 
power spectrum of the thick disk is well fit by a single power law PHi(f^) = AU~'^'^ which 
has the 3D slope. In this case we expect the break corresponding to the transition from 3D 
to 2D at f/ ~ 4k A. This break lies in the sample variance dominated region which explains 
why we do not detect it. 

We next incorporate the galaxy's radial profile using 
^/V2e ' 

where 9 = y/x'^ + y'^ is the radial coordinate in the plane of the disk. Here 6 coincides 
with the angle in the sky measured from the center of the galaxy. We have used 6q = 1' 
which corresponds to 128 mesh units . Note that the radial profile used in the simulation is 
exactly the same as the window function W{6) of the exponential model introduced in the 
previous section. The relative amplitude of ho and h{r) is a free parameter which decides 
the respective contributions from the smooth and the fluctuating HI components. Using the 



G{f) = exp 



expi-z'/zl), (14) 



ratio R = y {h'^)/hQ of the rms. value of h{r) to ho to quantify this, we have carried out 
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simulations for R = 0.35 and 0.7. We have also carried out simulations where the disk is 
tilted, and the normal to the disk make an angle of 60° to the line of sight. The simulated 
HI power spectra are shown in Figures [6] and [7] for the thick and thin disks respectively. 

We find that for a face-on disk with 6q = V, for both the thick and thin disks, the results of 
our simulations are in good agreement with our analytic estimate (Section 3) which predicts 
that the effect of the galaxy's radial profile is contained within a limited baseline range 
U < Um, which comes out to be f/ < Um ~ 3.5 9^^, for a = —2.5 (Figure [3]). The HI power 
spectrum is insensitive to the galaxy's radial profile at U > Um, and the shape of the power 
spectrum is the same whether we include the galaxy's radial profile or not. We find that the 
value of Um is not very sensitive to changes in R, the ratio of the fluctuating component to 
the smooth component. For a disk of the same size, the value of Um increases when the disk 
is tilted. The image of the disk tilted by 60° has an anisotropic window function W{6) with 
angular radius 1' and 0.5' along the major and minor axes respectively. For a tilted disk the 
value of Um is determined by the smaller of the two angular diameters, 0.5' in this case. 

In summary, our simulations demonstrate that for baselines U > Um the galaxy's radial 
profile has no effect on the power spectrum of HI fluctuations. In other words, the power 
spectrum would be unchanged if the same fluctuations were present in an uniform disk with 
no radial profile. Further, a break in the power spectrum at a f/ value greater than Um is 
indicative of a transition from 3D to 2D at the length-scale corresponding to the scale-height 
Zh = l/nU. 

Finally, we note that we have ignored velocity fluctuations and the galaxy's rotation 
throughout our simulations. This would rearrange the HI emission amongst the different 
frequency channels. Our simulation corresponds to the situation where there is just a single 
frequency channel whose width encompasses the entire HI emission, and the velocities have 
no effect in this situation. This is justified by the results presented later in this paper 
(Section 6) which show that, for the galaxies in our sample, the slope of the observed HI 
power spectrum does not depend on the width of the frequency channel. 



5 METHOD OF ANALYSIS 



The dwarf galaxy data used here is fro m Giant Metr e wave Radio Telescope (GMRT) ob 



ser vat ion s. The details can b e found in 
KK 230, begum et all feoosf ) for NGC 3741, 



Begum et al 



(2006) for UGC 4459, KDG 52 and 



Begum fc Chengaluil fl2003h for GR 8, Chen- 
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galur et al. 2008 (in preparation) for AND IV and lBegum fc Chengalui 



DuttaetaL 



(2004) for DDO 210. 
(120081 ) ■ for the spiral 



We use Very Large Array(VLA) archival data, discussed in 
galaxy NGC 628. 

n 

All the data are reduced in the usual way using standard tasks in classic AIPS Lj. For 
each galaxy, after calibration the frequency channels with HI emission were identified and a 
continuum image was made by combining the line free channels. The continuum was hence 
subtracted from the data in the u — v plane using the AIPS task UVSUB. 

The number of channels with HI emission (n) is different for each galaxy. To determine 
if the HI power spectrum changes with the width of the frequency channel, we combine N 
successive channels to obtain n/N channels and perform the power spectrum analysis for 
different values of in the range 1 < < n . The analysis is initially carried out for A^ = 1, 
and unless mentioned otherwise the results refer to this value. 

We measure Pm{U) by correlating every visibility Vp{U) with all other visibilities Vp{U + 
AU) within a disk | AU 9^-^ and then average over different U directions. Note that 
correlation of a visibility with itself is excluded. To increase the signal to noise ratio we 
further average the correlations in bins of U and over all frequency channels with HI emission. 

The measured visibility correlation estimator (eq. [T2|) is the convolution of the actual 
HI power spectrum with a window function. Here we do not attempt to estimate this win- 
dow function by fitting the observed HI image and then use this to deconvolve the power 
spectrum. Our approach is to estimate a baseline Um, such that for U > Um the effect of 
the convolution can be ignored. The measured visibility estimator is proportional to the 
HI power spectrum (eq. [13]) at baselines U > Um- We estimate the HI power spectrum us- 
ing only this range {U > Um) where the measured visibility correlation estimator may be 
directly interpreted as the HI power spectrum. 

Assuming an exponential window function, the value of Um depends on the parameters 
^0 and a. Note that we do not attempt to actually fit an exponential to the observed overall 
HI distribution and thereby determine 6q. The angular radius 6q of the galaxies in our sample 
is set to half the angular extent listed in Table [H The galaxies are typically ellipses and not 
circular disks as assumed in Section 3. Guided by our simulations, we have used the smaller 
of the two angular extents in our analysis. As Um oc 9q^, choosing the smaller value gives 
a conservative estimate of Um- The upper limit [/„ where we have a reliable estimate of the 



NRAO Astrophysical Image Processing System, a commonly used software for radio data processing. 
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HI power spectrum is determined by the requirement that the real part of PHi(f^) should be 
more than its imaginary part and also the noise in the line free channels. 

The procedure that we adopt to determine the best fit power law to the HI power 
spectrum is as follows. We first visually identify a baseline Um beyond which {U > Um) 
the visibility correlation estimator appears to be a power law in U . We then use the range 
Um < U < f/u to fit a power law PHi(f^) = A f/". The best fit a obtained by this fitting 
procedure is used to obtain a revised estimate for Um through Figure O The value of Um in 
this figure are for 6*0 = 1', we use the scaling Um oc O'^^ (as mentioned above, is taken 
to be the smaller value of the two angular extents for a particular galaxy) to obtain the Um 
value that is appropriate for the galaxy. The power law fitting is repeated with the revised 
estimate of Um- We iterate this procedure a few times till it converges. At each iteration, the 
best fit A and a were determined through a minimization. To test whether the impact 
of the window function is actually small, we have convolved the best fit power law with 
I Wv{U) p. The convolved power spectra were visually inspected to asses the deviations 
from the power law. The goodness of fit to the data was also estimated by calculating for 
the convolved power spectrum. The final fit is accepted only after ensuring that the effect 
of the convolution can actually be ignored. 

The values D/Um and D/Uu-, where D is the distance to the galaxy, give an estimate 
of the range of length-scales over which the power-law fit holds. In addition to the other 
parameters of the power law fit, these values are also shown in Table 2. 

6 RESULT AND CONCLUSION 

Figure [8] and [9] show the results of our analysis. The results are summarized in Table El 
We have detected the HI power spectrum of the galaxies DDO 210, NGC 628, NGC 3741, 
UGC 4459, GR 8 and AND IV. For these galaxies we find a range of baselines U where 
the real part of Phi{U) estimated from the channels with HI emission is larger than the 
imaginary part estimated from the same channels as well as the real part estimated from 
the line free channels. This is not true for KK 230 and KDG 52 (and also for AND IV with 
N = n), for these galaxies all three curves (Figure [9]) lie within the 1 — a errors-bars with 
very small offset and the interpretation is not straight-forward. For the subsequent analysis 
of all the galaxies we use only the real part of Phi{U) estimated from the channels with HI 
emission. 
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Galaxies 


DDO 210 


NGC 628 


NGC 3741 


UGC 4459 


GR 8 


AND IV 


KK 230 


KDG 52 


(1 a) N 


1 


1 


1 


1 


1 


1 


1 


1 


(2 a) velocity width (km s^-"^) 


1.65 


1.29 


1.65 


1.65 


1.65 


1.65 


1.65 


1.65 


(3 a) a 


-2.3 ± 0.6 


-1.6 ±0.2 


-2.2 ±0.4 


-1.8 ±0.6 


-1.1 ± 0.4 


-1.3 ±0.3 






(4 a) Um - Uu (kA) 


3.7 - 13.0 


1.0 - 10.0 


1.6 - 16.0 


3.1 - 22.0 


1.6 - 23.0 


0.6 - 6.7 






(5 a) D/Uu - D/Um (kpc) 


0.06 - 0.27 


0.8-8.0 


0.19 - 1.9 


0.18 - 1.16 


0.1 - 1.3 


0.56 - 6.2 






(6 a) x"/^ 


0.3 


0.2 


0.2 


0.7 


0.6 


0.4 







(1 b) N 


10 


64 


8 


12 


8 


16 


8 8 


(2 b) velocity width (km s~^) 


16.5 


82.56 


13.2 


19.8 


13.2 


26.4 


19.8 19.8 


(3 b) a 


-2.1 ± 0.6 


-1.5±0.2 


-2.5 ±0.4 


-1.7 ±0.4 


-0.7 ± 0.3 






(4 b) Urr^ - Uu (kA) 


3.7-13.0 


1.0 - 10.0 


1.6 - 16.0 


3.1 - 22.0 


1.6 - 23.0 






(6 a) D/Uu - D/Um (kpc) 


0.06 - 0.27 


0.8 - 8.0 


0.19 - 1.9 


0.18 - 1.16 


0.1 - 1.3 






(6 b) x^/>^ 


0.4 


0.2 


1.3 


0.3 


0.2 






(7) C2 


0± 1.2 


0±0.4 


0±0.8 


0± 1.2 


0±0.8 


0±0.6 





Table 2. The results for the 8 galaxies in our sample. Rows 1-6 gives (1) TV the number of channels averaged over, (2) the 
corresponding velocity width (3) a the best fit slope for the HI power spectrum, (4) U range over which the power law fit is 
valid, (5) length-scales over which the power law fit is valid and (6) the goodness of fit per degree of freedom. Row (7) 
Possible limits for the spectral slope of the velocity structure function. 



For all the galaxies we have tried to fit the observed Phi{U) with a power law Pm{U) = 
AU". We find that this provides a good fit with a reasonable x^/^ ^oi DDO 210, NGC 628, 
NGC 3741, UGC 4459, GR 8 and AND IV. The presence of a scale invariant, power law 
power spectrum indicates that turbulence is operational in the ISM of these galaxies. The 
range of length-scales for the power law fit differs from galaxy to galaxy and in total it covers 
lOOpc to 8Kpc. The details are summarized in Table O 

Both HI density fluctuations as well as spatial fluctuations in the velocity of the HI 
gas contribute to fluctuatio ns in the HI speciflc intensity. Considering a turbulent ISM, 



Lazarian fc Pogosyad (|2000| ) have shown that it is possible to disentangle these two contri- 



butions by studying the behavior of the HI power spectrum as the thickness of the frequency 
channel is varied. If the observed HI power spectrum is due to the gas velocities, the slope 
of the power spectrum is predicted to decrease with increasing thickness of the frequency 
channel. To test this we have repeated the power spectrum analysis increasing the channel 
thickness from iV = 1 to = n. In addition to A^ = 1, Figures [8] and [9l and Table [2] 
also show the results for N = n where the channel thickness spans the entire frequency 
range that has signiflcant HI emission. We do not flnd a significant change in the slope of 
the power spectrum for any of the galaxies. Since for all the galaxies the thickest channel is 
considerably larger than the velocity dispersion, we conclude that the HI power spectrum 
is purely due to density fiuctuations and not gas velocities. The fact that the slope does 
not change with channel thickness can be used to constrain the value of ^2, the slope of the 
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R (kpc) 



R (kpc) 



N = 1 



20 1 



N=1 



DD0210 




NGC628 




N = 10 




N = 64 



Real 

Imaginary 
Line Free 
Best fit 




0.5 1 



0.5 1 



T 0.5 



1 0.5 




10 20 



U (kA) 



U (kA) 



Figure 8. Power spectrum of the galaxies DDO 210, NGC 628, NGC 3741 and UGC 4459. The real and Imaginary parts of 
Phi {U) estimated after averaging N channels with HI emission, and the real part from N line- free channels are shown together 
for N = 1 (left panel) and N = n (right panel). The error-bars are for the real part from channels with HI emission. The best 
fit power law is shown in bold. In each case Um is marked with a bold-faced arrow and the fit is restricted to J7 > Um where 
the effect of the convolution with the window function can be ignored. 
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1 20 



N = 12 Real 

Imaginary 
Line Free 
Best fit 




U (kA) 



U (kX) 



N = 1 



AND IV 












N = S Real 










1 m a g 1 n a ry 










Line Free 






1 




Best fit 


























\ tI k 
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N = 1 



KK230 



N=8 Real 

Imaginary 
Line Free 



N = 1 



KDG52 



N=8 Real 

Imaginary 
Line Free 



U (kA) 



U (kA) 



Figure 9. Power spectrum of the galaxies GR 8, AND IV, KK 230 and M81DWA. The real and imaginary parts of Phi{U) 
estimated after averaging N channels with HI emission, and the real part from N line-free channels are shown together for 
AT = 1 (left panel) and N = n (right panel). The error-bars are for the real part from channels with HI emission. The best fit 
power law is shown in bold only where such a fit is possible. For those, Um is marked with a bold-faced arrow and the fit is 
restricted to f7 > Um where the effect of the convolution with the window function can be ignored. 
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velocity structure function (eg. iLazarian fc Pogosyanll2000l ). The C2 values are tabulated in 
Table El 

The galaxies in our sample have slope a ranging from —2.6 to —1.1. The two galaxies 

DDO 210 and NGC 3741 have slope 2.5, while the slope is 1.5 for NGC 628, 

UGC 4459 and AND IV, and —1.1 for GR 8. We h a ve pr oposed a possible explanation 



for this dichotomy in the values of a in 



Dutta et al. 



(I2OO8I ). This was based on the fact 



that DDO 210, where the power spectrum was measured across length-scales 100 — 500 pc, 
had a slope of —2.6 while NGC 628, a nearly face-on galaxy where the power spectrum 
was measured across length-scales 0.8 — 8kpc had a slope of —1.6. We have interpreted 
the former as three dimensional (3D) turbulence operational at small scales whereas the 
latter was interpreted as two dimensional (2D) turbulence in the plane of the galactic disk. 
For a nearly face-on disk galaxy we expect the transition from 2D to 3D turbulence to be 
seen at a length-scale corresponding to the scale height of the galaxy. Continuing with this 
interpretation implies that we have also measured 3D turbulence in NGC 3741, and 2D 
turbulence in UGC 4459, GR 8 and AND IV. 



The p ower spectrum analysi s can be used to dete r mine the scale- hei g. 



galaxies (lElmegreen et al. 



2001 



Padoan et al 



200 ih . 



Elmegreen et al. 



it of nearly face-on 



(I20011) have applied 



this method to estimate the scale height of LMC. In a recent study ( iDutta et al.l 120081 ) 
we have placed an upper limit of 800 pc for the scale-height of NGC 628. Our simulations 
(Section 4) show that we expect the slope of the HI power spectrum to change corresponding 
to a transition from 3D to 2D turbulence at the scale-height. The baseline U at which this 
break occurs can be used to estimate the scale-height D/tcU, where D is the distance to the 
galaxy. 

In addition to NGC 628, our sample contains 3 more nearly face-on galaxies with zhi < 
30°, namely DDO 210, UGC 4459 and GR 8. DDO 210 has 3D turbulence across the baselines 
2.0 — 10.0 kA. The absence of a break in the power spectrum places a lower limit of 160 pc 
on the scale-height. The galaxies NGC 628, UGC 4459 and GR 8 exhibit 2D turbulence for 
the entire U range (Table 2) in the measured HI power spectrum. This imposes the upper 
limit of 320 pc, 51 pc and 30 pc respectively on the scale- height . Note that the scale-height 
of 320 pc for NGC 628 differs from our earlier estimate (IDutta et al.l 120081 ) because of the 
extra factor of ir indicated by our simulations. 

The energy input fr om star formation is belie ved to be a major driving force for the 
turbulence in the ISM ( lElmegren fc Scald l2004al ). Hence, it is very interesting to check 
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a 




Figure 10. The slope a of the HI power spectrum plotted against Egpjj^, the SFR per unit area. The area has been determined 
from optical images and HI images in the left and right panels respectively. The galaxies with 3D and 2D turbulence are shown 
using empty and filled circles respectively. Note that the SFR for DDO 210, marked in the figure, is only an upper limit. The 
data for NGC 628, which is a spiral galaxy, is marked with an arrow. The straight lines show the linear correlation that we 
find (discussed in the text) between the slope of the HI power spectrum and the SFR per unit area. 

whether the slope of the power spectrum of intensity fluctuations in these galaxies has 
any correlation with the SFR. Any correlation, if present, will provide an insight into the 
processes driving tur bulence in the I SM. 

(120051 ) calculate the SFR per unit for 9 irregular 



In a recent study 



Willett et al. 



galaxies and investigate the correlation with the V-band and Ha power spectra. The length- 
scales they probe are 10 — 400 pc. In Ha they find that the power spectra becomes steeper 
as the SFR per unit area increases. However they do not find any correlation in the V-band. 
In this paper we probe length-scales 100 pc to 3.5 kpc. We test for a possible correlation 
between the slope of the power spectrum and the SFR per unit area and per unit HI mass for 
the 5 dwarf galaxies in our sample (Figure [TOl) . We report results using the area estimated 
from both, the HI disk and the optical disk. T he data for the SFR, angu lar extent, inc li natio n 
in HI and HI mass for these galaxies are from iBegum et al.l ( 120051 ) and iBegum et al.l (120081 ) . 
The SFR is determined from Ha emission, and the angular extent in HI is determined from 
the MomentO maps at a column density of 10^^ atoms cm~^. For the angu l ar extent and 



i nclina tion of the optical disk w e use the parameters from 



Young et al. 



(120001 ) 



Begum et al 



(12006! ) and lSharina et al.l (120081 ). The linear correlation coefficient is found to have values 0.35 
and 0.34 for the HI and optical disks respectively, indicating the absence of any correlation 
between the slope of the power spectrum and the SFR per unit area. We have similar result 
for SFR per unit HI mass for these 5 galaxies. However, we note a few points that should be 
taken into consideration when interpreting this conclusion. The Ha emission is not a good 



tracer of SFR in low mass galaxies due to stochastic star formation (IKarachentsev fc Kaisin 
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20071 : iLee et al.l 120071 ). Also, the length-scales across which the power spectrum has been 



estimated for the dwarf galaxies in our sample are substantially larger than the optical disk 
where star formation occurs. 

Turbulence can also possibly be related to different parameters of the galaxy. We test for 
correlations between the slope of the HI power spectrum and the following parameters: total 
HI mass, HI mass to light ratio, total dynamical mass, total baryonic mass, gas fraction, 



baryo n fraction. We use the estimates of these observable quantities from iBegum et al 



(120071 ). In all the cases the linear correlation coefficient is found to lie between —0.5 and 0.5, 
indicating the absence of correlation. 

In the above analysis we have considered galaxies with both 2D and 3D turbulence 
taken together, which in turn can suppress the correlations. Hence, we further investigated 
the correlations by considering the 4 galaxies in our sample (including NGC 628) with 2D 
turbulence. In this case the correlation coefficient comes out to be -0.70 and -0.74 for SFR 
per unit area of optical and HI disks respe c tively . This indicates a strong correlation and 



the result is similar to that of 



Willett et al 



(120051 ) . namely that the surface density of star 
formation rate is larger for galaxies with steeper HI power spectrum. These findings indicate 
a possible link between star formation and the nature of turbulence in the ISM. However, 
it would not be realistic to speculate on a possible cause and effect relation between these 
two. It is quite possible that the observed correlation is an outcome of extraneous factors 
which influence both. 

Finally, we note that the total number of galaxies in all of our analysis is rather small 
for a statistical conclusion. A larger galaxy sample is required for a better understanding of 
the generic features, if any, of turbulence in the ISM of faint dwarf galaxies. 
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